from subprocess import call

f = open("Read_DIV", 'w')

for j in [24,36,48]:
	for i in [20, 40, 60, 80, 100, 200, 400]:
		f.write('\n***threshold: %i maxLength: %i\n' % (i, j) )
		f.flush()
		call (["./mapSimulator", "human_g1k_v37.fasta", "human_g1k_v37.fasta.index.12", '-h', "testRead", "0", str(i), str(j)], stdout=f);
